Measurement of spatial qubits 



G Lima 1,2 , F A Torres-Ruiz 2 , Leonardo Neves 1 , A 
Delgado 2 , C Saavedra 2 and S Padua 1 

1 Departamento de Fi'sica, Universidade Federal de Minas Gerais, Caixa Postal 
702, Belo Horizonte, MG 30123-970, Brazil. 

2 Center for Quantum Optics and Quantum Information, Departamento de 
Fi'sica, Universidad de Concepcion, Casilla 160-C, Concepcion, Chile. 

E-mail: glimaOudec . cl 



Abstract. 

In this paper we study the state determination for composite systems of 
two spatial qubits. We show, theoretically, that one can use the technique of 
quantum tomography to reconstruct the density matrixes of these systems. This 
tomographic reconstruction is based on the free evolution of the spatial qubits 
and a postelection process. 
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1. Introduction 

The concept of quantum state plays a central role in Quantum Theory. The 
statistical distributions for the results of experiments carried out on a physical 
system can be completely predicted from its initial state. This has led to the 
development of techniques for the state determination. In the field of atomic physics, 
quantum endoscopy was used to determine the state of ions and atoms [TJ [2 [3] . In 
quantum optics, the Wigner function of multimodc fields could be measured using 
homodyne detection [HG3[5], and quantum tomographic reconstruction (QTR) was 
used for measuring the polarization state of the parametric down-converted photons 
[TJ, [SI [HI El HI]- Beside this, the raising of new technologic fields has brought more 
motivations for the study of these techniques of reconstruction. For example, in the 
field of Quantum Information, protocols like teleportation pj] and superdense coding 
|13| require initially knowing the quality of the quantum channel which is being used 
for implementing them. 

In spontaneous parametric down-conversion (SPDC) one photon from a pump 
beam incident onto a non-linear crystal originates, with small probability, two other 
photons usually called signal (s) and idler (i) or also twin photons [14]. Recently, we 
have demonstrated that by exploring the transverse momentum correlations of the 
twin photons it is possible to generate entangled states of two effective D-dimensional 
quantum systems [TS1 [TB] . Because the dimension of the Hilbcrt space of these photons 
is defined by the number of different ways available for their transmission through 
apertures where they are sent to, we call them spatial qudits. 

In this work we investigate the state determination of these quantum systems. 
We show, theoretically, that one can implement the process of QTR to obtain the 
density operators of two spatial qubits (D = 2) states. The theoretical description 
of our protocol is closely related to the description given in Ref. [5] for the quantum 
tomography of polarized two-qubit states. Even though we have considered just the 
case of qubits, the description used here can be generalized to spatial qudits states. 

2. The state of two spatial qubits 

The angular spectrum of the pump beam is transferred to the two-photon 
state generated by SPDC, when the monochromatic, paraxial and thin crystal 
approximations are assumed [T71HU 19, 20, 21]. The spectral matching is implemented 
by inserting narrow bandwidth interference filters in front of the detectors [17[ I21j . 
The state of twin photons when they are transmitted through generic apertures, under 
above approximations, can be written in one dimension as [151 117] 



where | lqj } is the Fock state for one photon in mode j = s,i with transverse 
wavevector qj. T{q s , qi) is the two-photon amplitude and it is given by 



The longitudinal distance za defines the plane where the apertures are placed, and 
k p is the wave number of the pump beam used to generate the twin photons. The 




(1) 




x W(±(x s + Xi); z A ) exp (-i(q s x s + qiXi)). 



(2) 
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function W(x; za) is the pump field distribution at the plane of the apertures (z = z A ) 
and at the transverse position x. Aj(xj) is the transmission function of the aperture 
in mode j. 

In the case where these apertures are two identical double slits, the two-photon 
state in Eq. |T|) simplifies to [15j [16] 

\V) = W+, + \ + ) a \+).+W + ,-\ + ) a \-) i 

+ W-,+ \ - )„ | +> < +W_,_| -) s \-) v (3) 

where W+- = W-.+ = 1 W(0;z A )e ik " d2 / 8zA and W±,± = jW(±/2; z A ). d is the 
separation between two consecutive slits and a is each slit's half width. 7 is a 
normalization constant. The | + )j (or | — }j) state is a single-photon state defined, up 
to a global phase factor, by 

\±).=^J, dq 3 e^ ld sine ( qj a)\ \q 3 ), (4) 

and represents the photon in mode j transmitted by the upper (lower) slit of its 
respective double slit. These states form an orthonormal basis in the two-dimensional 
Hilbert space of each photon and are used to define the twin photons' logical spatial 
qubit states. Therefore, it is clear that the state | \& } represents a system composed 
of two spatial qubits. 

A typical setup used to generate the spatial qubits is outlined in Fig. [TJ A pump 
beam incident upon a non-linear crystal generates SPDC the twin photons which are 
correlated in their transverse momenta. These photons are then sent through identical 
double slits (A; and A s ) placed at the z^-plane and, after being transmitted by these 
apertures, the twin photons will be in the state of Eq. (J3]). 




Figure 1. Schematic diagram of the experimental setup used for generating 
spatial qubits. BBO is a non-linear crystal used to generate the twin photons. 
A s and Ai are the double-slits placed at signal and idler propagation paths, 
respectively. 



3. Detecting the spatial qubits in distinct bases 

For doing the quantum tomographic reconstruction of two-qubit states codified 
in polarization it is necessary to consider the detection of these qubits in some 
polarization states. The measurement bases that are usually considered in the 
experiments are [Sj: {\H}, \V)} (horizontal and vertical polarizations), {|D),|A)} 
(linearly polarized light at 45° and —45° with respect to the horizontal polarization, 
respectively), and {\L), \R)} (left and right circularly polarized lights, respectively). 
The measurement (the projection) in one of the states of these bases is represented by 
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a projector operator, and an adequate choice among the available projections allows 
the reconstruction of the composite state. 

Now we show how the spatial qubits can be measured in distinct bases of 
{|±) s ,|±)i}. After being transmitted by the double slit, the twin photons freely 
evolve in space. Because of diffraction, it spreads out faster in the x-direction so that 
it passes from the discrete states | ± } , to states continuously distributed along this 
transverse direction. In a certain z-plane situated far behind the double slits' plane 
(z = za), each particular superposition of the basis states |±) will have a specific 
transverse probability distribution [22] , To recover the discrete nature of the logical 
states in this z-plane, we need to implement an adequate postelection process for the 
twin photons. As we shall show in the following lines, this can be properly done by 
allocating in this plane one double slit for the transmission of each mode j. We refer to 
these new apertures as detection double slits. The transverse position of their slits are 
Xj k where the subscript j denotes the mode at which the double slit was placed (again 
j = s, i) and k = 0, 1 stands for the lower and upper slit of this aperture, respectively. 
The slits' width is 2b. The superscript /j, is used to define a set of measurements. For 
a certain value of [i wc have two transverse positions available, Xj and Xj 1; for the 
slits of the double slit. It can be seen as the mean position of the double slit's center 
and, as we shall see, it is the index which defines the spatial measurement basis which 
is being considered. The state of the twin photons that crossed these two detection 
double slits can be written as 

J2 ^i/(< fc )>i/(< fc o>. (5) 

k,k'=0,l 

where the states \f(x^ k ) ) are the states of the postselected photons which crossed the 
slit in the transverse position xf t of the double slit placed in mode j. These states 
are given by 



( x 1kb \ 
x sine I qj + qjb I | lqj ) , (6) 

where a = (z — ZA)/k p and one can see that the logical states used to describe the 
twin photons are once again discrete. The calculation of the state | $} is shown in 
|Appendix A| However, what is really interesting to note now is that there is a relation 
between the coefficients of this state and the coefficients of the original state of the 
twin photons given by Eq. That is 

Bk.k> = ^2 r u{x^ k )r v (x^ k ,)W u ^ v , (7) 



with 



j,k T d) a 
r ±( x ^ k ) = sinc ( V - 2q ] ex P 



V 



4a 



(8) 



Therefore, one can think of the state of the twin photons transmitted through the 
detection double slits [Eq. ([5])] as the original state of these photons [Eq. J3])] rewritten 
in the new basis {\f(x^ k )),\f(x^ k ,))} (again k, k' = and 1). As wc mentioned 
before, the spatial bases are designed in terms of \x. For distinct values of \l one 
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has distinct transverse positions x^ and x^ 1 for the slits of the detection double 
slit in mode j. The states \f(x^ ) ) and \ f(x^ x ) ) are orthornomal, i.c, they form an 
orthornomal basis for the Hilbert space of the single photon after the detection double 
slit. We use the superscripts /i, //, fi" , ... to denote distinct bases. The post-selection of 
the twin photons in one of the basis states of {\f(Xg k ) ), \f(xf k ,) )} is then represented 
by the following projector: 

P <Kl = I /«*), /«*)></«*), /«*) I- (9) 
4. Reconstructing the qubits state 

As is discussed in Ref [8] , the density operator of two qubits is specified by sixteen real 
parameters. Thus, for determining their state, one needs to construct sixteen equations 
that are linearly independent and that are functions of measurable quantities. From 
now on we refer to these equations as MLI equations (measurable and linearly 
independent equations). In this section we show how this can be done for the spatial 
qubits state of Eq. ([3]). To do this we first remind the reader of some basic results of the 
two photons interferometry theory. The photodetection of a s-photon (at r,t) followed 
by the detection of a i-photon (at r',t') is described by the fourth-order correlation 
function which is defined as [2] 

(r,r<, t,t') = ($\E { -\r, t)E^ (r\ t')E^ (r' , t')E\ +) (r, t)| * >, (10) 

where | $ } is the two photon state and E^ are the positive and negative parts of the 
electric field operator. The coincidence rate of these photons (at r,t,r',t') is, therefore, 
given by 

Cr.rw =£G (2) (r,r',i,0, 

where £ is a constant related to the total number of photon pairs sent to the detection 
apparatus and to its efficiency. 

4-1. The protocol 

We can now explain the protocol for doing the reconstruction of the two spatial qubits 
state given by Eq. (|3|). To start with this, let us first assume the following general 
formula for its density operator p: 

= ^ Pi,m-i'^i'\l s ;m l )(l' s ;m' i \, (12) 

l,m,V ,m' 

where the indices l,m,l',m' assume the values + (or — ) when denoting the upper 
(lower) slits of the double slits Aj's. For example, the matrix element, is the 

probability of detecting the twin photon signal at the slit "+" and the idler photon 
at the slit "-" . 

To construct the first set of MLI equations, one should first measure the 
coincidence rates with two detectors {Di and D s ) placed just behind the slits of the 
apertures A4 and A s as is shown in Fig[2]^a)[j]. Four MLI equations can be constructed 

{ For practical purposes, it is worth mentioning that this measurement can also be done in the plane 
of image formation of these apertures I23| . 



Measurement of spatial qubits 



6 



by permuting the detector positions behind the slits of these two double slits. The 
coincidences rates behind these slits, according to Eq. (fTTj) . can be written as 

C Um = r)NTr(pP l>m ), (13) 

where N is the number of the photon pairs transmitted through the double slits and r\ 
the coincidence counts efficiency for the detectors used. P; j7ra is the projector operator 
given by P; jjn = \ l si m i )(l s ,m i \. If we work out this expression we will see that 

C;,m = V N Pl,m;l,m: ( 14 ) 

and, therefore it is clear that, after doing the normalization of the coincidence counts 
seen in this part of the experiment, one will be able to determine the diagonal elements 
of the spatial qubits density operator 




Figure 2. Schematic diagram of the spatial measurements needed for doing the 
quantum tomography of two spatial qubits. In (a) the coincidence counts are 
recorded with the detectors (D s and Di) fixed behind the slits of the double slits 
A s and A;, respectively. In (b) they are fixed behind the detection double slits 

whose slit transverse positions are: Q = and x\j ^ = In (c) one detector 
records the idler photons transmitted by the slits of A; and the other, the photons 
transmitted by the signal detection double slit. 

To obtain four more MLI equations one should now introduce the detection double 
slits in the z-plane, with their slits at the transverse positions x^ = and x 1 ? 1 = ^ 
(again, a = (z — ZA)/k p ), and record the coincidence rates with the detectors Dj's 

§ Here is considered the usually accepted fair sampling assumption where the ensemble observed in 
the experiment (r) X N twin photons) is considered to be a fair sample of the total ensemble of p (TV 
twin photons transmitted by the double slits). 
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behind these slits [See Fig.[^b)]. The new coincidences rates are now written in terms 
of the projectors of Eq. © and are given by 

f2b\ 2 

x ,x 1 \ _L J x o ' x i 

f2b\ 2 

<V*- = T vNTr(pP , ), 
x i ' x o \ Li J i ' 

C .f ,.f = ( r J ^(pP^, xf 0> (15) 

where L is the transverse size of the idler and signal diffraction patterns (caused 
by the apertures Aj's) at the z-plancjjj]. Therefore, it is clear that the factor ^ 
is a normalization constant which accounts for the effect of one photon diffraction, 
which causes a reduction of ^ in the detector's efficiency that is being used to record 
this photon at the z-plane. From these measurements one will be able to determine 
the imaginary and real parts of some of the matrix elements of p. In Table [TJ we 
give a summary showing which are the components of p that are determined from 
the measurements considered in this protocol. There is also a comparison with the 
measurements done in the QTR of two polarized qubits [8] . In | Appendix B| we give 
the expressions which relate the elements of p and all the coincidence rates considered 
in this protocol. 

Let us now consider the detection of the signal photons behind its detection double 
slit at the z-plane and the detection of the idler photons behind the slits of the double 
slit Ai (at the z^-plane) as it is schematically shown in Fig. E^c). The four possible 
coincidence rates obtained, again by permuting the detectors' positions behind the 
slits of these double slits, are now written as 

c x ?', m = T r]NTr (f )P ^,J> ^ 

with again m = + (m = — ) for the detection of the idler photon at the upper (lower) 
slit of Ai. The projector P x ^ m is defined by 

k ' 

p x», m = \f(x^ k ),m){f(x^ k ),m\. (17) 

These other four equations are MLI equations and they allow the determination of 
new elements of p as is shown in Table [TJ 

The last four MLI equations are obtained when the signal photons are detected 
behind the slits of A s and the idler photons behind the slits xf = and ^fi = 

|| This is, of course, an approximation where we assume that the total length of the double slit 
diffraction pattern, L, is defined by the transverse length of the three principal peaks of diffraction. 
However, we can say that Eq. il 51 is a good approximation for the values of the coincidence rates 
considered since most of the intensity of the light scattered by a double slit lies in this L region. 
K Here it should be clear that we are always considering that the MLI equations constructed are 
linearly independent in the system of sixteen equations that we want to construct. The reader can 
check this by calculating Eq. J7J, explicitly, for the values x^ Q = and x^ x = ^ used. 
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Reconstructing p 



Measur. 


Pol. proj. 


1 


\HH) (HH\ 


2 


\HV) (HV\ 


3 


\VH) (VH\ 


A 


|T/1/\ /T/T/ 
| V V J \ V V \ 


5 


\HD) (HD\ 


6 


\HL) (HL\ 


7 


\VD) (VD\ 


8 


\VL) (VL\ 


9 


\DH) (DH\ 


10 


\LH) (LH\ 


11 


\DD) (DD\ 


12 


\DL) (DL\ 


13 


\DV) (DV\ 


14 


\LV) (LV\ 


15 


\LD) (LD\ 


16 


\LL) (LL\ 




Table 1 



Coefficients 



Spatial proj. 



Coefficients 



Ph,h-,h,h 

PH,V;H,V 

PV.H-V.H 

PV,V;V,V 
R {PH,H;HV) 
1(PH,H;H,v) 

R (pv,V;V,H) 
^(PV,V;V,h) 

R (ph,h-,v,h) 

I(PH,H;V,h) 

R (PH,H;V,v) j R {PH,V;V,h) 
1 {PH,H;V,v) , I (PH,V:V,h) 

R (iOV,V;H,v) 

I(pv,y ; if,y) 
I {PH,H:V,V) , I (ph,V:V,h) 
R (ph,h-,v,v) , R (ph,V:V,h) 



P+.+ 
p+.- 
P-.+ 

P- _ 

P- x" 

p „< 
p M' , 



/>+,-;+ - 
P-,-,-,- 
%+,+;+,-) 

I(P+,+;-,+) 

"(/=>+,-;-,+) 

%-,-;+,-) 
i (p +1+; _,_) > i(p + _ 1 _ + ) 

i(p + , +; _,_),R(p +) _ ; _,+) 



R(p +I+ i-,-), 



tomography of two spatial qubits. There is also a comparison with the 
measurements used in the QTR of polarized qubits. \H), \ V), \D) and \L) 
are the kets representing polarized qubits with horizontal, vertical, diagonal and 
left-circular polarization, respectively. The indices R and I stand for the real and 
imaginary parts, respectively. 



at the transverse z-plane. The projectors which represent these four coincidence 
measurements are given by 

P»^ = |i,/(< fc )><f,/(< fc )|, (18) 

and the coincidence rates are 

2b 

Ct, x » = -r,NTr(pP lixS ), 

c i, x f = T vNTr{pP ^f ] ' (19) 

where I = + (I = — ) holds for the detection of the signal photon at the upper (lower) 
slit of A s . The elements of p that can be determined from these last measurements 
are also shown in Table [1] 

So, after performing the sixteen measurements described above, one will have 
reconstructed the density operator p of the two spatial qubits state generated in the 
experiment considered. It is worth mentioning that the protocol presented here can 
also be applied to any type of two spatial qubits state, including the mixed states 
recently reported in the literature [24] . 

Here, we give a brief description of an experimental setup that can be used for 
implementing this reconstruction scheme. This description is based on already known 
parameters from our experiments. First, to generate spatial qubits (qudits) we use 
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two identical double (multi) slits aligned in the direction of the signal and idler beams, 
which are positioned at a distance z a = 200 mm from the Type II BBO crystal. The 
slits's width is 0.1 mm with a separation of 0.25 mm. The photons transmitted 
through the double-slits are detected in coincidence. In front of the detectors there 
are interference filters, with very small bandwidths (with a typical FWHM of 1 nm) 
and centered at twice the pump beam wavelength, which will select the frequency of 
the generated signal and idler photons [171 12 lj . For reconstruction purposes identical 
single slits at the detection plane can be used at position z — z a = 600 mm. At 
this distance, the fourth order interference pattern has a width of the order of 8 mm. 
The detection slit has a width of 0.1 mm. This slit can be easily installed in front 
of single photon detector modules. If a Kripton laser, with an average power of 
200 mW is used as a pump beam, a maximum number of 200 coincidence counts in 
500 s are expected for (b) configuration in Fig. 2. For this experimental setup we have 
Xq = and Xi = 0.496 mm. These measurement positions are within the fourth order 
interference pattern region and they can be well resolved. In case of the measurements 
of type (a) and (c) 200 coincidence counts are expected for measurement times of 20 s 
and 100 s, respectively. Furthermore, there are not requirement for numerical aperture 
of the detectors. 

4-2. Maximum likelihood estimation 

The technique of QTR is based on a linear inversion of the measured data as was shown 
above for the two spatial qubits states. Thus, it is dependent on any experimental 
errors that may occur while recording the data. They can appear as a consequence of 
experimental noise or misalignment and the reconstructed state is, in general, only a 
reasonable approximation of the real quantum state. The density matrices obtained 
may have properties that are not fully compatible with a quantum state. 

To generate only possible density matrixes there is an alternative that has been 
used during the state determination. It is the numerical technique called maximum 
likelihood estimation [23]. However, even though it generates only possible density 
matrixes, it has the drawback of enhancing the uncertainty in the state estimation. 

Once the two qubit density operator is reconstructed with the technique proposed 
here, one can apply the matrix obtained to this numerical optimization. In Rcf [5] a 
practical approach for doing this is extensively discussed and it can also be applied 
for two spatial qubits systems. 

5. Conclusion 

In conclusion, we have theoretically shown that one can use the technique of quantum 
tomography to reconstruct states composed of two spatial qubits. Even though we 
have considered the state determination only for this special case, we believe that 
it can be generalized and performed in a similar way for spatial states of more 
dimensions. The importance of this work comes from the possibility of using the 
spatial qubits for quantum communications protocols where the ability to characterize 
may be necessary. The main advantage that we can envisage at this stage is the 
possibility of using transverse correlations of the down-converted photons for encoding 
quantum information in qudits instead of qubits. Besides, even in the case of qubits 
we see some other important advantages. For instance, when entangled polarization 
qubits are propagated through optical fibers they suffer decoherence because of the 
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depolarization effect of the fibers. The birrcfrigcncc of the fiber gradually destroys 
the entanglement of the entangled polarized qubits. The entanglement of the time-bin 
entangled qubits also suffers of severe decoherencc in optical fiber due to the chromatic 
dispersion effect. However, when the spatial qubits are used for sending information 
thought optical fibers, both depolarization and chromatic dispersion does not affect, 
directly, their quantum correlations and thus the entanglement of the spatial qubits 
seems to be more robust against the decoherence effects of optical fibers. We have 
already started this line of research in our group to demonstrate this idea. 
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Appendix A. Calculating the state | <3> ) 

As we showed in Section [3J the state of the down-converted photons when they are 
transmitted through double slits, placed at the z^-plane, is given by Eq. ([3]). Now, 
if we consider the free space propagation of the twin photons from these apertures to 
another transverse z-plane, we will obtain a new form for their state [23j : 

l*)„ « [dq s /"^(%,g s )|lft)|l? s ), (A.l) 
where the new two-photon amplitude Q (qi,q s ) is 

Q {q s ,qi) = ?-e~ laq h- taq s sine (%a) sine (q s a) (w + ,+e- t{q > +q ° )d + W/_._e lte+<Zs)d ) 

! sine (qia) sine (q s a) (V+^e" 1 ^-^ + W-,+e l ^- q ° )d ^ , 

(A.2) 



a 2 

TT 



(z — Za) 

again, a = - k ■ 

Now we consider the determination of the twin photon state after they are 
transmitted through the detection double slits considered in our protocol and which 
are placed now at this transverse z-planc. The transverse position for the slits of these 
double slits are Xj and Xji- To calculate this state we first assume a general form 
for it 



|$) cx J dq s J dqiD (q s , q t ) \lq s ) \lq t ) . (A.3) 
The two-photon amplitude T> (q s , qi) can be written in terms of the convolution 

V(q s , qi ) = fdq'JdqiGiCq^Tsiqi - q^T^ - (A.4) 



where T^'s are the Fourier transform of the transmission function of the detection 
double slits and are given by 

T( qj ) = 2b\e iq '<o S inc ( qj b) + e iq ' x ^ sine ( qj bj\ , (A.5) 
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^^-^W Sine 
sine ((q' s - q s )b) ^ W u , v e« +V ^) x £ e <(?:-?.)<* e *M-*> 



u,v— ± 



fc.fc'=0,l 



By using the following solution 

I± ( x » kiqj )=[ e-'^'e^V^^^smc (gjo) sine ((^ - % » dgj 



se 3 J - fc e 4q sine 



and by denning the coefficients 



r±(rf <k ) = sine 



2a 



cxp z- 



2q 



2a 



(A.7) 



4a 



(A.8) 



one can write the function T> (q s , qi) as 

T>(q s ,qi) = ^ E Wu > vru ( x s,k) r v( x i.k>) e ~ lqXi ' k ' sinc 



fe,fe'=0,l tt,-u=± 



x e I<?;E = fc sinc 



2a 



2a 



qb 



qb 



(A.9) 



Now if we replace Eq. (|A.9j) into Eq. (|A.3|) and define the new slit states 













2a 


x sinc ^ 



M |i<&>, 



(A.10) 



we get that the state of the twin photons which were transmitted by the detection 
double slits is 



$)cx Yl E ^«« fc )^(< fc ') /«,,)) (A.n) 

fe,fe'=0,l u,u=± 

= ^ Sfc,*|/(<j)|/(<*)), 

fc,k'=0,l 



with 



u,v=zL 



(A.12) 
(A.13) 



(A.6) 
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Appendix B. The coefficients of p 



Here we show explicitly the expressions for determining all the coefficients of the 
density operator, p, by following the measurement scheme given in Section 14.11 [Sec 
Fig. [5]). The first type of measurement is performed with the detectors Di and D s 
right behind the slits of the double slits A; and A s , respectively (See Fig. [2ja)). As 
it was discussed above, the diagonal coefficients are given by 



V N Cl ' m ' 



(B.l) 



with the normalization that ijN = C+.+ + C+,_ + C_,+ + C_,_. We remind that Ci, m 
is the coincidence rate recorded with the detector D s behind the slit I of its double 
slit and with the detector Di behind the slit m of Aj. 

Now, we consider the detection of the signal photon behind its detection double 
slit at z-plane [Sec Fig. [UJc)] while the other detector remains behind the slits of the 
aperture A s (See Section |4~T1 for more details). In this case the coincidences rates are 
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C x ^ m = —r]NTr(pP x ^ m ), 

C < \ m = TVNTr(pP x ^J, 
and by considering the values of =0 and 1 = %j we get 
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where \ — L/(2b). 

Reversing the role of signal and idler detectors, by considering the coincidence 
rates of Eq. (fH)l) in our protocol, we obtain 

P+,+;+- = —Tt ( X fc+.xfl - iC, »>) + ° + ' + t ° + '~ ( _1 + J = Pa 
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The last elements of p to be calculated are the anti-diagonal coefficients. To 
determine them we need to consider the measurements which involve the propagation 
of both photons, i.e, the coincidences rates of Eq. (fT"5)) . The anti-diagonal elements 
are given in terms of these coincidences by 
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